Simulating fermion production in 1 + 1 dimensional QED 
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We investigate fermion-anti-fermion production in 1 + 1 dimensional QED using real-time lattice 
techniques. In this non-perturbative approach the full quantum dynamics of fermions is included 
while the gauge field dynamics can be accurately represented by classical-statistical simulations 
for relevant field strengths. We compute the non-equilibrium time evolution of gauge invariant 
correlation functions implementing 'low-cost' Wilson fermions. Introducing a lattice generalization 
of the Dirac-Heisenberg-Wigner function, we recover the Schwinger formula in 1 + 1 dimensions in 
the limit of a static background field. We discuss the decay of the field due to the backreaction of 
the created fermion-anti-fermion pairs and apply the approach to strongly inhomogeneous gauge 
fields. The latter allows us to discuss the striking phenomenon of a linear rising potential building 
up between produced fermion bunches after the initial electric pulse ceased. 

PACS numbers: ll.lO.Kk, 11.15.Ha, 12.20.Ds 



I. INTRODUCTION 

It has already been pointed out in the early days of 
quantum physics that the vacuum of quantum electrody- 
namics (QED) becomes unstable against the formation 
of many-body states in the presence of strong external 
electromagnetic fields, manifesting itself as the creation 
of electron-positron pairs by the Schwinger mechanism 
[ll-Q. Nevertheless, this fundamental quantum effect has 
not been experimentally observed so far as it has not 
been possible to generate the required electromagnetic 
field strengths in a laboratory. However, due to the rapid 
development of laser technology during the last decades 
an experimental verication of electron-positron pair pro- 
duction in the focus of high-intensity laser pulses comes 
into reach. 

Vacuum pair production in an applied uniform elec- 
tric field of strength Eq may be viewed as a quantum 
process in which virtual electron-positron dipoles can be 
separated to become real pairs once they gain the bind- 
ing energy of 2mc^. However, there will be strong spa- 
tial and temporal inhomogeneities of the electromagnetic 
field in realistic situations as envisaged in upcoming high- 
intensity laser experiments. The theoretical description 
of such a non-perturbative phenomenon in quantum field 
theory out of equilibrium is a demanding task and very 
little is known so far for realistic scenarios. Most current 
approaches assume the electromagnetic field as being an 
external one with a one-dimensional inhomogeneity, so 
that the problem of particle production can be mapped 
onto a one-dimensional quantum mechanical scattering 
problem [1,[H. This approach neglects, in particular, the 
backreaction of the created fermion-anti-fermion pairs on 
the electromagnetic field. This is closely related to kinetic 
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descriptions in terms of a momentum dependent distri- 
bution function of pairs in collisionlcss (Vlasov) approx- 
imations [6l-[lll|. For multi-dimensional inhomogeneities, 
more advanced approaches such as semi-classical approx- 
imations of the vacuum effective action [l2| or the Dirac- 
Heisenberg-Wigner phase space formulation [l^, [T^] have 
been appHed. However, to describe strongly inhomoge- 
neous field configurations including the full backreaction 
of the produced particles remains a theoretical challenge. 
In view of the potential experimental applications it is 
crucial to devise new theoretical methods which can deal 
with this situation. 

In this work we propose to use real-time lattice gauge 
theory techniques [l5| to compute fermion-anti-fermion 
pair production in QED. In this non-perturbative ap- 
proach the full quantum dynamics of fermions is included 
while the gauge field dynamics can be accurately rep- 
resented by classical-statistical simulations for relevant 
field strengths. As the inclusion of dynamical fermions 
can become numerically very expensive, the real-time 
evolution of fermions is taken into account by means of 
a low-cost fermion algorithm [l6j . 

As an example and in order to compare with estab- 
Hshed continuum results, we apply these techniques to 
QED in 1 + 1 dimensions - the massive Schwinger model 
[itI [TsI . Introducing a lattice generalization of the Dirac- 
Heisenberg-Wigner function, we show that the simula- 
tions accurately reproduce the results described by the 
Schwinger formula in the limit of a static background 
field. We discuss the decay of the field due to the backre- 
action of the created fermion-anti-fermion pairs and ap- 
ply the approach to strongly inhomogeneous gauge fields. 
For these fields we compute for the first time the back- 
reaction of the created pairs on the gauge fields. Most 
strikingly, we find that a self-consistent electric field be- 
tween the produced fermion and the anti-fermion bunch 
builds up for times exceeding the initial pulse duration. 
The two bunches consisting of particles and anti-particles 
act as a capacitor, creating a homogeneous electric field 
between them, which can be represented in terms of a 
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linear rising potential. 

From the point of view of theoretical developments, 
it is important to note that very similar questions are 
addressed in physics of the early universe or in collision 
experiments of heavy nuclei. Non-equilibrium particle 
creation from large coherent fields has been extensively 
studied in the context of inflaton dynamics [l9l - [2H and 
non-Abelian gauge field theory [22|, The produc- 

tion of fermion-anti-fermion pairs has been mainly in- 
vestigated based on semi-classical descriptions using the 
Dirac eq uation coupled to time-dependent background 
fields j24| - |27| . Going beyond these approximations, stud- 
ies based on two-particle irreducible (2PI) effective ac- 
tion techniques [H, [2^ showed that quantum effects can 
dramatically affect the far-from-equilibrium production 
of fermion pairs. Recently, pair production from scalar 
inflaton decay [l^ as well as baryogenesis [s^, HH has 
been studied using similar lattice field theory methods 
as employed in this work. In particular, the results of 
Refs. [m provide a proof of principle that real-time 
lattice simulations with Dirac fermions are indeed fea- 
sible in 3 -I- 1 dimensions. They also have been tested 
[29| against calculations based on real-time 2PI effective 
action techniques [H, H^l in their range of applicability. 

This paper is organized in the following way: In Sec.HIl 
we briefly review the low-cost fermion algorithm and de- 
rive the real-time lattice equations of motion for the mas- 
sive Schwinger model. Additionally, we construct a lat- 
tice generalization of the Dirac- Heisenberg-Wigner func- 
tion which is subsequently used as a read-out tool for 
fermionic distributions. In Sec.|lll]we first apply this for- 
malism to a static electric background field and compare 
to the Schwinger formula. We then discuss the decay of 
the background field due to the backreaction of the cre- 
ated fermion-anti-fermion pairs. As a second example, 
we investigate the pair creation process in the presence 
of a space- and time-dependent electric field. In Sec. HVl 
we conclude and give an outlook. 



II. REAL-TIME LATTICE GAUGE THEORY 

A. Continuum formulation 

We consider QED in 1 + 1 dimensions, which is defined 
in the continuum by the action 

5 = Jd^x (^^[i^D^ - m]i; - ^-^'^'-^m-) , (1) 

with the covariant derivative ~ dfj_ + leA^ ensuring 
gauge invariance of the action under local U{\) transfor- 
mations 

^^Ve*^^ , A^^A^-d^,K. (2) 

Here = 0, 1 as space-time is only two-dimensional with 
x^ = t and x"^ = x. The field strength tensor F^'^ = 



d^A'^ — d'^A^ possesses only one non-trivial component 
which is regarded as the electric field: 

= -J^^ = E{x,t) . (3) 

We will frequently consider temporal axial gauge with 
Aq{x, t) = Q and simply denote the spatial component of 
the vector potential as A{x,t)} One observes that the 
electric field E[x,t) is the canonical momentum conju- 
gate to A[x, t). 

The Dirac algebra is composed of two Dirac gamma 
matrices only: 

= 2g^"' with (7^)^ = ^'^-^^'^^ , (4) 

with g^" = diag(l, —1). This algebra may be represented 
in terms of the first two Pauli matrices 7*^ = cti and 
7^ = ~i<72. Moreover, the chirality matrix 

{7^^,75} = with (75)^ = 75 , (75)' - 1 , (5) 

can be defined in terms of the third Pauli matrix 75 = 
(T3. As a consequence, the spinors ^ and ■0 are two- 
component field operators, obeying the equal-time anti- 
commutation relation: 

mx,t),^{y,t)} = j°S{x-y) . (6) 

1. Time evolution equations 

In general, in the classical-statistical theory observ- 
ables are calculated as ensemble averages of solutions of 
Maxwell's equation 

d,T^''ix,t)^{rix,t)) (7) 

starting from different canonical field variables at initial 
time to, here Ata{x) ~ A{x,tQ) and Eti,{x) ~ E{x,to). 
The values for the canonical field variables at initial time 
arc distributed according to a normalized phase-space 
density functional WjAtg, Etg], such that an observable 
(O) is given by [13, Hi: 

(O) = J DAt,DEt, W[At,,Et,]0,^[At,,Et,] . (8) 

Here 0,Mt„,Et„] = ^ DAO[A]5{A - A,AAt,,Et,:\), 
where ^ci[-4tQ, is the solution of the classical field 
equation ([7]) with initial conditions = At^ and E^i — 
Eta at initial time t^. Ensemble averages at initial time 
are taken to correspond to the respective quantum expec- 
tation values for the gauge fields. The gauge field dynam- 
ics in the classical-statistical approximation is accurately 
described in the presence of sufficiently high occupation 



We note that this incomplete gauge choice leaves a residual gauge 
invariance under time-independent gauge transformations. 
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numbers or fields, which is in general the case for the rel- 
evant field strengths for pair production. It breaks down 
once the typical gauge field occupancies become of order 
unity. For an introductory review see Ref. js^ . 

The subsequent time evolution then follows from ^ 
with 



(9) 



where the expectation value is taken with respect to the 
initial state of the spinor field. We will restrict ourselves 
to the Dirac vacuum within the current investigation. 
The equations of motion for the spinors read: 



iD*^%l){x,t)^^ = —mtp(x,t) 



(10a) 
(10b) 



Since the fermions appear only quadratically in the action 
(H]) these equations are exact for given classical gauge 
field configuration. 

Equivalently, the above equations can be conveniently 
expressed in terms of the equal-time statistical propaga- 
tor 

Fix,y;t) = ^{[4,ix,t),^iy,t)]) , (11) 
which yields the closed system of equations: 



ij^'D,,^F{x,y;t) 
,yF(a;,y;t)7^ 



iD: 



d^F>'''{x,t) 



mF{x,y;t) , (12a) 
-mF{x,y;t) , (12b) 
-eTi[j''F{x,x;t)] . (12c) 



We note that the gauge field dynamics in 1 -I- 1 dimensions 
is special since it is governed by the fermionic backreac- 
tion only. Therefore, we do not consider sampling over 
initial gauge field configurations in this work. 



2. Initial conditions and low-cost fermions 



We have to solve the Cauchy problem (|T2|) in order to 
calculate fermion-anti-fermion pair production. Accord- 
ingly, we need to provide an initial value for the statis- 
tical propagator at to = 0. To this end, we consider an 
asymptotic Dirac vacuum - corresponding to zero parti- 
cle number and vanishing gauge field - and employ the 
framework of canonical quantization: 

J 27r 

(13) 

with ui = \/ rn? -|- p^ and anti-commuting creation and 
annihilation operators 

{a(p),at(p')} = {6(p),&^(p')} = 27r5(p-p') , (14) 



whereas all other anti-commutators vanish. An explicit 
representation of the eigenspinors is given by 



u{p) 
v{p) 



1 +p 
^2uj{uj+p) V 

1 fuj +p 



(15a) 
(15b) 



^2uj{oj+p) V 

fulfilling the orthogonality relations: 

{p)u{p) = 1 = v'^{p)v{p) , v) (j))v{—p) = . (16) 

Because of the fact that the asymptotic Dirac vacuum 
is homogeneous in space and time, we obtain the initial 
value: 

^,.,,,,^/|.*-'>=^. (17) 

The solution of the time evolution equation of the sta- 
tistical propagator ([6]) may be based on a mode function 
expansion [15i] . This treatment can be well suited for low 
dimensional systems but becomes computationally too 
expensive in higher dimensions. In view of later applica- 
tions of our approach to 3 -I- 1 dimensional systems, we 
perform a stochastic integration of an equivalent set of 
equations going by the name of low-cost fermions p^ . 
To this end, we introduce ensembles of classical stochas- 
tic spinors, termed male ipM^Xji) and female ippix^t), 
instead of spinor field operators. Given these c-number 
spinors, we define: 

Fsto{x,y;t) = {iljMix,t)ipFiy,t)) = {il;F{x,t)ipM{y,t)) , 

(18) 

where here (...) is understood as an ensemble average. 
The requirement 



Fstoix,y]t) = F{x,y;t) 



(19) 



is met provided that the stochastic spinors ipg^x^t), with 
the gender index g = {M, F} both satisfy the Dirac equa- 
tion (|10ap and Fsto{x,y;t) takes the initial value (fTT]). 
This second requirement is achieved by initializing the 
stochastic spinors according to 



M = j ^e^^^l=\u{p)^{p)±v{^p)^{p)] , (20) 



^g{x 



with complex random variables ^(p) and 7?(p). Note that 
the male and female spinors only differ by the sign of the 
antiparticle component. In order to reproduce the initial 
value (|17[) . the random variables are sampled according 
to 

{me ip)) - Hp)v {P)) = i2n)S{p -p'), (21) 

whereas all other correlators vanish. 

In an actual simulation employing low-cost fermions, 
the closed system (fT2l) is solved in the form: 



ij^'Df,ipg{x,t) = rmljg{x,t) , (22a) 
d^T'"'{x,t) = -eTr[7''^^,to(^,^;t)] • (22b) 
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The stochastic spinors ilig{x,t) are evolved in time in- 
dependently and the ensemble average (...) appearing in 
the definition (fT5|) is approximated by an average over 
a sufliciently large number A'sto of pairs of male/female 
spinors. While the computational cost of the mode func- 
tion approach scales with the volume of the phase space, 
i.e. N'^'^ in d dimensions, the resource requirements of 
'low-cost' fermions are proportional to just N^Nsto- 



B. Lattice formulation 



We solve the equations of motion ([22|) on a 1 -I- 1 
mcnsional space-time lattice. For the spatial sublatt: 
we define: 



1 di- 
;ice, 



A 



e {0,...,7V,-i} 



(23) 



with the spatial lattice spacing as and the total number 
of spatial lattice sites Ng- A point on the space-time 
lattice is then denoted by x = with the temporal 

lattice spacing at such that t = atj- We employ periodic 
boundary conditions in the compactificd spatial direc- 
tion whereas no periodicity assumptions apply for the 
non-compact temporal direction. The lattice action gov- 
erning the dynamics is then composed of a pure gauge 
part as well as part describing the fermions, including 
their interaction with the gauge field. 



1. Pure gauge part 

In order to put the gauge fields on the lattice, we use 
the compact formulation of a gauge theory with C/(l) 
symmetry. The parallel transporter Ufj,{x.) is associated 
with the link from a lattice point x to a neighboring 
point X -I- /i in the direction of the space-time lattice axis 
/x = 0,l: 



„ieaj,.A,,(x) 



(24) 



The link variable obeys C/^(x) = J7~^(x) and we use the 
definition C/_^(x) = f/*(x — /i). The continuum gauge 
transformation ([5]) translates to 

?/'(x) ^ 0(x)V'(x) , (25a) 
C/p(x) ^ r!(x)C/^(x)ri*(x + /i) , (25b) 

with O e U{1). Given the gauge-dependent link variable, 
we define the gauge-invariant plaquctte variable: 



[/^,(x) = [/^(x)[/,(x + A)[/;(x + i>)[/:(x) 



(26) 



Disregarding higher order terms in the lattice spacings 



a^, we find: 



f/^.(x) 



(27) 



Accordingly, the pure gauge part of the action can be 
written as 

[U] = 5^ Re [1 - f/oi (x)] . (28) 

X 

Moreover, the electric field ([3]) is given by 

1 



i?(x) 



■Im[C/oi(x)] 



(29) 



2. Dirac and interaction part 

Using a symmetric finite difference approximation 
for the first derivatives, the naive discretization of the 
fermionic part is given by 



4"^[^,^",C/] =a,a,^V^(x) 

X 

U + A) - ^^^(X)V^(X - A) 



2an 



■ mi/'(x) 
(30) 



where the gender index is omitted for simplicity. This 
expression is gauge-invariant under lattice gauge trans- 
formations ((25|) . however, it also gives rise to unphysi- 
cal states. This fermion doubling problem is well-known 
from Euclidean lattice field theory. However, unlike to 
Euclidean time one distinguishes between temporal and 
spatial doubler modes in the real-time formulation [isj . 

The spatial doublers, corresponding to high- 
momentum excitations showing a low-energy dispersion 
relation, are conveniently suppressed by adding a higher 
derivative term to the action: 



(31) 



which vanishes in the continuum limit as — > 0. This 
gauge-invariant Wilson term in space ensures that only 
low-momentum excitations show a low-energy dispersion 
relation. In the lattice implementation, this corresponds 
to adding one more term to the action: 



S^^^ [ijj, ijj, U] = atas ^ ■i/'(x) 

X 

"C/i(x)V'(x + i) - 2?A(x) + [/_i(x)V'(x - i) 



2a., 



(32) 



We do not include a temporal Wilson term as this would 
turn the Dirac equation into a second order differential 
equation in time. The temporal doublers are avoided 
provided that we initialize only the physical mode and 
choose the temporal lattice spacing to be much smaller 
than the spatial lattice spacing at as [HI, [l^, . 



5 



3. Low-cost Wilson fermions 

The construction of the stochastic spinor ensemble on 
the space-time lattice follows the same lines as in the 
continuum outlined in Sec. Ill Al For spatial momenta p 
we define the conjugate lattice: 



A 



Lp 
2tt 



2 



(33) 



with the spatial volume L = asNs- Accordingly, the 
discrete Fourier transformation is given by 



qeA 



(34a) 



^5]exp(-^')V'(x), (34b) 



the complex random variables ^{q) and r]{q) are sampled 
according to 



This is most easily done by assuming 



(40) 



= X{q)e^^^^^ , rj{q)=Y{q)e 



ie(q) 



(41) 



and choosing the amplitudes X{q) and Y{q) to be 
Gaussian distributed whereas the phases 4){q) and 9{q) 
are chosen to be uniformly distributed on the interval 
[-7r,7r). 

4-. Lattice equations of motion 



with the notation q = (q^j)- The stochastic spinors are 
then again initialized at to ~ at jo = 0; 



1 



^s(qo) = -^[u{q)C{q) ± vi-q)r]{q)] 



(35) 



To simplify simulations afterwards, we use the gauge 
freedom and employ the lattice equivalent of the tempo- 
ral axial gauge: Uq (x) = 1 for the equations of motion. 
Stationarity of the lattice action 



with qo = {q,jo) and eigenspinors: 

1 



S[4',^,U]=Sg[U]+ 5f [^, V, U] + >[^,^,U] (42) 



u{q) 
v{q) 



UJ + q 

^2Cj{Cj + q) y 

1 / LU + q 



y^2CbiCb + q) V 



(36a) 
(36b) 



Here the mass term is modified due to the spatial Wilson 
term 



with respect to the temporal link Uo{x.) results in the 
discrctizcd version of the Gauss law: 



E{x) - E{x - 1) = ^V(x)7V(x + 0) -t- c.c. (43) 



2 . 2 f '^q 

m — m -\ sm — 

a, V Ns 



and we define the lattice quantities 



~ 1 . /27rg\ 
q = — sm 



(37) 



(38) 



In order to reproduce the correct initial value for the 
statistical propagator on the lattice. 



f 27ri(/i — h)q\ m — qj^ 



iV, 



2uj 



(39) 



This equation is a constraint which is fulfilled during the 
time evolution for the considered initial conditions. 

The stationary condition of the action with respect to 
the spatial link C/i(x). on the other hand, results in the 
equation of motion: 



£;(x)-£;(x - 6) = 



eat 



- i]f/i(x)V7(x + i) + c.c. (44) 



Finally, the stationarity condition of the action with 
respect to the Dirac field '0(x) gives: 



V'(x + 6) = V(x - 6) - 2tat (m + — ) 7V(x) - — - i]f/i(x)V'(x + i) - + z][/_i(x)V'(x - i)) 

\ as J as ^ 



(45) 



The set of equations ([^51) - (|^5|) is the lattice version of (1^^ in temporal axial gauge including a spatial Wilson 



6 



term. 

In order to solve the Cauchy problem, we have to pro- 
vide the following initial values at to — at jo ~ 0: 



E{^o-0) , [/i(xo) , ^(xo-6) 



^(xo) 



with xo = {I, jo) for all I E A. Most notably, we have to 
choose initial values for the spinors at Jq — 1 and jo, which 
is a consequence of the chosen leapfrog algorithm. To be 
able to initialize them we assume a free field evolution at 
initial times. 

The algorithm, which is a variant of the one introduced 
in [36j . can then be summarized in the following way: 

1. Electric field evolution: Given i?(x — 0), Ui{x.) and 
■0(x) we evolve the electric field to E{x.) according 
to (1441). 



2. Dirac field evolution: Given ip{:x. — 0), V'(x) and 
t/i(x) we evolve the Dirac field to ^(x + O) accord- 
ing to (|45|) . 



3. Temporal plaquette: We evaluate the temporal pla- 
quette C/oi(x) according to (|27|) : 



C/oi(x) 



(46) 



4. Spatial link evolution: The link variable C/i(x -I- 0) 
is calculated from the temporal plaquette Uoi{x.) in 
temporal axial gauge according to p6)) : 

C/i(x + 6) = ;7oi(x)C/i(x) . (47) 

5. Reiterate the steps 1-4. 

C. Gauge-invariant correlation functions 

In order to compare our simulation results with typi- 
cal discussions u sing the Dirac- Heisenberg-Wigner phase- 
space approach IsTIT^. ISTt - lilj . we define suitable gauge 
invariant two-point correlation functions on the lattice. 



1. Continuum Wigner function 

Starting from the continuum expression for the sta- 
tistical propagator a gauge-invariant generalization 
may be defined as: 

F{xi,X2',t) = cxp [ ie / dxA{x,t) \ F{xi,X2;t) . 



(48) 

The Wilson line factor ensures gauge invariance under 
local C/(l) transformations. The Fourier transformation 
with respect to the relative coordinate defines the Wigner 
function: 

W{x,p, t) = ~ I dye-'Pyp {x + y/2, x - y/2; t) , (49) 



with X = {xi -\- a:2)/2 and y = Xi — X2- The Wilson line 
factor in (|48|) is not unique, however, a physical sensible 
interpretation of p as kinetic momentum forces the inte- 
gration path to be chosen along the straight line. Equiv- 
alently to (j49p . we may also write: 



yV{x,p,t) 



dze^'P'-''-'^Fiz,2x- z;t)+j.c. (50) 



with the abbreviation: 

D + 'j.c. = D + j°D^j° 



(51) 



As the Wigner function is in the Dirac algebra and fulfills 
yyt = 7°yV'7°, one can decompose it in terms of its Dirac 
bilinears: 



yy = i [s + i75p + 7°TVo - 7^ 



(52) 



where all its irreducible components can be chosen to be 
real. Regarding the Dirac vacuum, which is described by 
the statistical propagator p7|) . the only non- vanishing 
components are given by: 

Svac(a:,p,t) = , ^viic{x,p,t) = . (53) 

(jj (jj 

In terms of these components the total charge Q and 
the total energy £ can be expressed as phase-space inte- 
grals: 



Q = e / dTwo{x,p,t) 



(54a) 

£ = y dr[ms{x,p,t) + p-v{x,p,t)] + — J dxE^{x,t) , 

(54b) 

with the phase-space volume element dT = dxdp/{2TT). 
The integrands e{x,p,t) = [ms{x,p,t) + pv{x,p,t)] and 
g{x,p,t) = vo(x,p, t) are regarded as energy pseudo- 
distribution and charge pseudo-distribution, respectively. 
We may define further quantities such as the particle 
number pseudo-distributions: 



n^(x,p,t) 



e{x,p, t) - evac(cc,p, t) ± uJwo(x,p, t) 
2uj 



(55) 

which may be associated to the density of particles and 
anti-particles, respectively. Of course, in the interacting 
quantum theory the interpretation of these phase-space 
pseudo-distributions, collectively denoted as m{x,p,t), 
has to be taken with care. We emphasize that our ap- 
proach is not based on these quantities and we use them 
only for read-out and comparison with literature results. 
We will frequently consider also the partially integrated 
position space and momentum space marginal distribu- 
tions: 



mx{x,t) 



dp 
2^ 



m{x,p, f) 



mAvA = \dxmi^x,p,t) 



(56a) 
(56b) 
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or the fully integrated quantities: 

m(<) = J drm{x,p,t) , (57) 
instead of the pseudo-distributions m(x,p,t). 

2. Lattice Wigner function 

In order to adjust the above continuum treatment to 
the lattice, we have to account for the periodicity of the 
spatial lattice properly. Our approach is an extension 
of previous work on the discrete Wigner function in the 
context of signal processing . 

We first define the gauge invariant generalization of 
the lattice statistical propagator according to: 



(58) 



where U{li,l2',j) is the lattice analogue of the Wilson 
line factor along the straight line path. However, since 
the straight line path between two lattice points is not 
unique due the periodicity of the lattice, we choose to de- 
fine it such that properties of the above standard contin- 
uum interpretation apply. It turns out that this requires 
taking the shortest path between two lattice points. Ac- 
cordingly, for Al — li — I2 > we employ: 



Al < 



Al>^ 
2 



U=l[ (7i(x) , 

l = l2 

U= U C/r(x) X n C^i*(x) 



l=li 1=0 

On the other hand, for < we use: 

/2-1 



Al > 



2 



Al<-f 



u=Y[ [/r(x) , 

Ns~l h-l 

Yl Uii^) X Y[ C/i(x) 



(59a) 
(59b) 

(60a) 
(60b) 



More precisely, we utilize the following Wigner lattices: 



Aw = 
Aw = 



- e{0,...,27V,-l} 

a,. 



TT 



(61a) 
(61b) 



which have the same extent as the original ones A and A, 
however, each with twice as many grid points. We then 
define the lattice Wigner function according to 

W(/,g,.7)^-^e"V^= 



^Y^e'^-^k,/N^Fik,[l~kU/,j)+l.c. 

(62) 



feeA 



with I £ Aw and q € Aw- We account for the periodic- 
ity of the lattice by taking the module operation in the 
second argument of the statistical propagator: 

[I ^k]N,^ {I- k) mod Ns . (63) 

This definition is such that we reproduce the above con- 
tinuum expressions for the marginal distributions, as 
shown in Appcndix|Al Moreover, the lattice Wigner func- 
tion dM]) again fulfills W'^ = 7°W7° so that the decom- 
position in terms of its Dirac bilinears ((52|) is possible. 

In complete analogy to the continuum, we may then 
again define various pseudo distributions: 

(64a) 
(64b) 

(64c) 



g{l,q,t) = e^o{l,q,t) , 
e{l,q,t) = [rhs{l,q,t) +qv{l,q,t)], 
±n „ -^\ _ e(',9,0-evac(Z,g,t)±(Ins/o(/,(7,t) 



n {l,q,t) 



2Cj 



corresponding to charge, energy and particle/anti- 
particle number, respectively. Given these pseudo- 
distributions m{l, q, t), the marginal distributions are de- 
fined via 



mxihj) = W7 ^ rn{l,q,j) 



2L 
2 



mviqj) = ^ rn{l,q,j) , 



(65a) 
(65b) 



ieAv 



whereas the fully integrated quantities are given by 

Here one should note the summation order in the last 
expression: The sum over I G Aw yields the marginal 
distribution m-p{q,t) which is non- vanishing for even q 
only. Accordingly, the subsequent sum is just taken over 
qG A. 

III. PAIR PRODUCTION SIMULATIONS 

We now come to the results which are based on the 
lattice approach presented in the previous section. As 
a first example, we consider a static electric background 
field, disregarding the backreaction of created fermion- 
anti-fermion pairs. This configuration can be solved an- 
alytically such that we can compare our lattice simu- 
lations with well established continuum results. Sub- 
sequently, we also include the backreaction of created 
fcrmion-anti-fermion pairs and discuss the decay of the 
gauge field which shuts pair production off after a char- 
acteristic time. 

As a second example, we investigate the pair creation 
process in the presence of a space- and time-dependent 
electric field. Neglecting backreaction in a first step, 
we can compare to and complement previous investiga- 
tions based on the continuum Dirac-Hcisenberg- Wigner 
approach [iO, |4^ . Subsequently, we solve the full lattice 
evolution and compare. 
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A. Spatially homogeneous gauge field 

We consider a static electric background field E{x, t) = 
Eq in temporal axial gauge Aq = 0, represented by the 
vector potential 



A{t) = Eot 



(67) 



Within the compact lattice formulation, this corresponds 
to a trivial temporal link L'o(x) = 1 and the spatial link 



C/i(x) = e"'''"''''^"^ 



(68) 



disregarding higher order terms in the lattice spacing. 
Moreover, we introduce the dimensionless field strength 
parameter 



Eo 
Ec ' 



(69) 



with the critical Schwinger field strength Ec ~ rr? je. For 
all subsequent numerical results we employ ejm = 0.3. 
In Appendix [B] wc briefiy review some analytic results, 
which are used for comparison in the following. 



1. Particle production without backreaction 

In this section we disregard the backreaction of created 
fermion-anti-fermion pairs on the electric field. This cor- 
responds to neglecting the fermionic contributions in the 
gauge field equation of motion pi)) . Starting with the 
vacuum initial conditions for the fermions, this amounts 
to evolving the fermion equation ([45|) with a sudden 
switching-on of the electric field at initial time. 

In Fig. [T]we show the time evolution of the total num- 
ber of produced particles, n~^{t), for various values of the 
dimensionless field strength parameter e. Most notably, 
we observe two different regimes: At early times there is 
a transient oscillatory behavior superimposed which can 



be attributed to the sudden switching-on of the electric 
field. For e = 1 we estimate this oscillation to be ex- 
ponentially damped with a characteristic rate 7 ~ 1/m, 
leading to a purely linear growth to very good accuracy 
after times of a few 7"^. 

The slope of the linear rise of n+ (t) strongly depends 
on the value of e. In order to extract its functional depen- 
dence, we perform a linear fit. For this we measure the 
change in the total number of particles An^ which are 
produced during the time interval T = 10/m for times 
large compared to 7^^. In Fig. [2] we compare the slope 
for different values of e with the analytical result from 
Appendix |B] 



An+ e 



(70) 



We emphasize that for this analytical result the initial 
time is sent to the remote past such that it cannot re- 
produce the transient oscillatory regime. However, both 
the simulation and the analytical result should accurately 
agree for large enough times. 

The lattice results are shown in Fig. [5] for different 
spatial lattice spacings keeping the volume L = a^iVs 
constant, thus increasing Ns accordingly. One clearly 
observes that the simulation and the analytical result 
([70)1 fall nearly on top of each other for small enough 
as, indicating that we are close to the continuum limit 
in that case. As a matter of fact, we find that temporal 
discretization errors are quite negligible for at < 0^/20. 
This corroborates that the real-time lattice simulation is 
in fact capable of reproducing the analytic results in the 
continuum limit to very good accuracy. 

In Fig. [3] we show the normalized particle number 
marginal distribution n^{p,t)/L, corresponding to the 
momentum spectrum of created particles (|B2|) . and com- 
pare it to the continuum value f{p). In comparison to 
the integrated particle number shown above, the spec- 
trum is not smooth but shows fiuctuations due to the 
sampling of low-cost fermions. As a matter of fact, these 
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e= 1.0 




6 = 0.9 
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6 = 0.6 
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FIG. 1: Time evolution of the total particle number n'^{t) for 
dilferent values of e. The parameters are A'sto = 10^, at — 
0.00125/m, as = 0.025/m, A^^ = 1024 such that L = 25.6/m. 



FIG. 2: Comparison of the analytical results (|70|) with the 
numerical fit for L = 25.6/m and dilferent lattice spacings Os- 
The remaining parameters are A'^sto = 10'^ and at — as/20. 
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«p(p,t)/L 




p/m 



FIG. 3: Comparison of the continuum expression /(p) (solid 
line) with the normalized particle number marginal distribu- 
tion n^{p,t)/L (dashed line) for e — 1 at t = 50/m. The 
parameters are A^sto = 10'', at — 0.00125/m, as = 0.025/m, 
Ns = 1024 such that L = 25.6/m. 




FIG. 4: Time evolution of the total particle number n^(t) 
with (solid line) and without (dashed line) backreaction for 
an initial value e — 1. The parameters are A^sto ~ 10**, at = 
0.0025/m, as = 0.075/m, iV^ = 256 such that L = 19.2/m. 



fluctuations can be systematically reduced by taking A^sto 
larger. We find that it suffices to take the number N^to 
of the order of 10'^ in order to accurately calculate in- 
tegrated quantities such as n+(t). This is in contrast 
to the momentum spectrum rip (p, t) where the number 
A'sto needed to be at least of the order of 10"' to suppress 
the statistical fluctuations sufficiently and obtain sensible 
results. In contrast to the one-dimensional case consid- 
ered here, the convergence is expected to be even better 
for three space dimensions where self-averaging plays a 
major role [29| . 

The established interpretation of f{p) is such that elec- 
tric field energy is taken and transformed into virtual 
fermion-anti-fermion pairs, showing up as the distinc- 
tive peak around momenta p = 0. If the applied field 
strength Eq is large enough, i.e. of the order of Ec, these 
charged excitations can be separated over the Compton 
wavelength and become real fermion-anti-fermion pairs. 
These real particles are then further accelerated in the 
background electric field and achieve higher and higher 
momenta up to p — > oo. 

We observe good agreement of simulation and ana- 
lytical results regarding the virtual fermion-anti-fermion 
peak around p = as well as the overall magnitude of 
n^{p,t)/ L. However, we observe a qualitatively different 
behavior for large momenta. This is due to the fact that 
the analytic result assumes an electric field which has ex- 
isted for all times such that all momenta up to p — > oo 
are already occupied whereas we solve an initial value 
problem on the lattice. Accordingly, we observe a tran- 
sient effect corresponding to the peak at high momenta 
propagating to higher and higher momenta during the 
time evolution. 



2. Particle production with backreaction 

We now include the backreaction of created fermion- 
anti-fermion pairs on the electric field. As a consequence, 
particle creation comes with a simultaneous decrease of 
the electric field due to energy conservation. This en- 
ergy transfer from the gauge sector to the fermion sector 
finally results in a decay of the electric field. 

In Fig. [4] we compare the time evolution of the total 
particle number n^{t) for simulations with and without 
backreaction. We have already seen in the previous sec- 
tion that the particle number grows eventually linearly 
if we disregard the backreaction of created fermion-anti- 
fermion pairs on the electric field. However, this changes 
drastically if we include the backreaction mechanism: 
Following the transient regime at early times, the pair 
production rate immediately slows down once the pair 
creation process kicks in and the electric field is weak- 
ened. Eventually, this is getting to a point where the 
fermion-anti-fermion production process effectively stops 
and n+ (t) levels off. This process happens on rather short 
time scales of the order of At ^ 25/m. 

To see the long-time behavior, in Fig. \5\ we show the 
particle number n~^{t) and the electric field E{t) for times 
up to 800/to. Most notably, we observe the occurrence 
of plasma oscillations in accordance with previous inves- 
tigations 0: Starting from to = 0, the magnitude of the 
electric field decreases due to the creation of fermion- 
anti-fermion pairs. Due to the backreaction mechanism, 
an internal electric field builds up so that the field eventu- 
ally changes sign and grows until a first local minimum is 
achieved. The electric field then increases again, changes 
sign, reaches a local maximum and so forth. The oscilla- 
tion frequency fl increases with the number of produced 
fermions, in accordance with the expected parametric de- 
pendence. 

The behavior of the particle number n"^ (t) follows from 
the oscillatory behavior of the electric field: Particle ere- 
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40 200 400 600 800 



t-m 

FIG. 5: Time evolution of the electric field E{t) (dashed line) 
and the total particle number n'^{t) (solid line) for an initial 
value e = 1 for a much longer time period. The parameters 
are as in Fig. |4] 



ation effectively terminates when the magnitude of the 
field strength drops below ~ 0.5Ec, corresponding to the 
approximate plateaus in n~^(t). However, at those in- 
stants of time at which the electric field reaches local ex- 
trema, fermions are created again. Due to the fact that 
the envelope of the electric field decreases with time, the 
particle number n~^{t) assumes the shape of a staircase 
with decreasing step height. 

We emphasize that the classicality condition (AA) ^ 1 
[s^ l is well fulfilled also after the backreaction effectively 
terminates the pair production: For an electric field am- 
plitude E with characteristic oscillation frequency the 
classicality condition reads E'^/fl^ ^1. In our case 
E ~ Ec/2 ~ 771^ /2e during these times such that with 

~ ttto/SO for the employed coupling e/m 0.3 we 
have E^/n^ ~ 700. 

Moreover, in Fig. [6] we demonstrate that the energy 
transfer from the gauge sector to the fermion sector is in 



agreement with energy conservation. 

Finally, in Fig. [7] we compare the normalized parti- 
cle number marginal distributions (p, t) /L for simula- 
tions with and without backreaction prior to the onset of 
plasma oscillations. We observe two major modifications 
if we include the backreaction: 

First, the high-momentum peak is shifted to lower mo- 
menta. This is due to the fact that acceleration in an 
electric field is proportional to its field strength. Ac- 
cordingly, particles are less accelerated and achieve lower 
momenta if the electric field is decreasing gradually. 

Second, the overall magnitude of n^[p,t)/L declines 
in the low-momentum regime. Again, this can be at- 
tributed to the decay of the electric field as the decrease 
of the field strength is accompanied by a drop in the pair 
production rate. Consequently, this gradual decrease of 
the pair production rate shows itself as a decreasing am- 
plitude of n^{p,t)/ L. 

This simple picture changes rather drastically at late 
times because of the occurrence of plasma oscillations. In 
Fig. [8] we show the normalized particle number marginal 
distribution n^{p, t)/L at different times. Due to the fact 
that the electric field changes its sign again and again, 
the fermions are accelerated back and forth in momentum 
space over and over again. The shaking of the fermions 
by the electric field has several implications: 

In contrast to the wedge-shaped spectrum at early 
times, this results in a peaked n^{p, t)/L at late times. It 
has to be emphasized, however, that this peaked distribu- 
tion still oscillates around p = in accordance with the 
electric field. Moreover, owing to the ongoing creation 
of fermion-anti-fermion pairs at times when the electric 
field reaches its local extrema, the overall magnitude of 
n^{p,t)/ L increases as well. 




200 



400 
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p/m 



FIG. 6: Energy transfer between the fermionic sector (solid 
line) and the gauge sector (dashed line) for an initial value 
e — 1. The dotted line shows the total energy, with the 
fermion vacuum contribution being subtracted. The parame- 
ters are as in Fig. (4] 



FIG. 7: Normalized particle number marginal distribution 
rip (p, t) /L with (solid line) and without (dashed line) backre- 
action for an initial value e = 1 at t = 50/m. The parameters 
are A^sto = 10^ at = 0.00125/m, a, = 0.025/m, iV, = 1024 
such that L = 25.6/m. 
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B. Space- and time-dependent field 

As a further example we consider an inhomogeneous 
electric background field which is localized in space and 
time: 



(71) 



E{x,t) = _BoSCch (a;i)cxp ( 



where w and A determine the duration and spatial ex- 
tent of the pulse, respectively. Studies based on the con- 
tinuum Dirac-Heisenberg-Wigner function only recently 
started to address such inhomogeneous configurations, 
disregarding the backreaction of created fermion-anti- 
fermion pairs [1^ H^j ■ Here we are for the first time able 
to take this fermionic backreaction into account using our 
lattice techniques. This will allow us to discuss the strik- 
ing phenomenon of a linear rising potential building up 
between produced fermion bunches for times exceeding 
the pulse duration. 



1. Particle production without backreaction 

In a first step, we solve the problem without taking into 
account backreaction. Consequently, we do not evolve 
the electric field according to (|44|) as it does not fulfill 
Maxwell's equation. We rather force the electric field to 
be given according to (|7ip at every space-time point and 
investigate the fermion-anti-fermion production in this 
given background field. 

In Fig.[5]wc show the position-space marginal distribu- 
tions n^{x,t) for three different times, with the electric 
field parameters e = 1, w = 0.1m and A = 5/m. One ob- 
serves two qualitatively different regimes, corresponding 
to early times ('creation regime') and late times ('prop- 
agation regime'). 

The fermion-anti-fermion pair creation process takes 
place at early times, when charged excitations are created 

4(p,t)/L 




t=50/m 



p/m 



FIG. 8: Normalized particle number marginal distribution 
rip (p, t) /L for simulations with backreaction at different times 
for an initial value e = 1. The parameters are A^sto = 10^, at ~ 
0.0025/m, as = 0.075/m, Ns = 512 such that L = 38.4/m. 



in a space region where the electric field acts. The cre- 
ation process also comes with a polarization effect, sep- 
arating positive from negative charges. It has to be em- 
phasized, however, that n'^{x, t) and n^{x, t) still overlap 
at these early times. 

This changes in the propagation regime: Owing to the 
acceleration by the electric field, one bunch of excita- 
tions with positive charge propagates into the positive 
x-direction whereas another bunch of excitations with 
negative charge propagates into the opposite direction. 
Asymptotically, these bunches can be identified with par- 
ticles and antiparticles, respectively. 

In Fig.[Tn]wc show the total number of created particles 
ri+(t) for t — > oo as a function of the spatial extent A. The 
result without backreaction corresponds to the dashed 
line. One clearly observes the termination of the fermion- 
anti-fermion creation process for small values of A: The 
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FIG. 9: Position space marginal distributions n^{x,t) (solid 
line) and n J (p, t) (dashed line) for e = 1 at different times 
t = (top), t = 0.6/a; (middle) and t = 1.2/aj (bottom). 
The parameters are Nsto ~ 10^, at = 0.01/m, as — 0.22/m, 
Ns = 256 such that L = 56.32/m. 
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pair creation process terminates if the work done by the 
electric field over its spatial extent is too small to provide 
the rest mass energy of the fermion-anti-fermion pair. 
This observation is in perfect agreement with previous 
studies [H, H m ■ For large values of A we find a 
linear growth of the particle number, which reflects the 
scaling of the available electric field energy that grows 
with A. 



2. Particle production with backreaction 

We now consider the numerical solution of the full lat- 
tice problem including backreaction. The solid line in 
Fig. [TO] shows the full result for the total number of cre- 
ated particles — )■ cx)) as a function of the spatial 
extent A. In accordance with the previous discussion we 
find that the backreaction mechanism tends to decrease 
the number of created particles since the electric field is 
weakened by the pair-production. For large enough spa- 
tial extent of the pulse, such that the pair-production is 
significant enough for backreaction to become important, 
this eventually decreases the slope of the curve growing 
linearly with A for large spatial extent. 

The dashed curve in the upper part of Fig. [TT] shows 
the position-space marginal distribution n^{x,t) at time 
t = Q/uj. At this time the duration of the initial electric 
field pulse and the corresponding pair creation regime 
is long over. The electric field parameters are e — 1, 
w = 0.2m and A = 5/m. 

The acceleration by the electric field leads to one bunch 
of excitations with positive charge propagating into the 
positive x-dircction whereas another bunch of excitations 
with negative charge is propagating into the opposite di- 
rection. Most strikingly, we find that a self-consistent 
electric field E{x,t) between the two fermion bunches 




A-m 



FIG. 10: Total number of created particles n {t oo) for 
e = 1 as function of the spatial extent A of the pulse. Com- 
pared are the full result (solid line) and the result without 
backreaction (dashed line). The parameters are A'sto ~ , 
at = 0.00125/m, Os = 0.22/m, iV^ = 512 such that L = 
112.64/m. 





0.4 




U.J 




0.2 


a 


0.1 









-0.1 




0.04 




U.U J 




0.02 


t)/m 




>< 


0.01 


+i>< 


c 







-0.01 





FIG. 11: Self-consistent electric field E(x,t) (solid line) and 
position space marginal distributions ri^ {x,t) (dashed line) at 
t — Q/u for e = 1, ui = 0.2m and A = 5/m. The arrows indi- 
cate the propagation direction of the particle and anti-particle 
bunch, respectively. The parameters are A'^sto ~ 10^, at = 
0.0075/m, as = 0.22/m, iV^ = 512 such that L = 112.64/m. 



builds up in the absence of any external field ([71]) at these 
times. The two bunches consisting of particles and anti- 
particles act as a capacitor [i^, creating a homogeneous 
electric field between them whereas there is no field out- 
side them. This electric field is shown in the lower part of 
Fig.[TTJ Owing to the description of the fermionic degrees 
of freedom in terms of low-cost fcrmions, we observe some 
small fiuctuations in the electric field on top of this ho- 
mogeneous field. Again, these fiuctuations decrease with 
increasing A^sto- 

The homogeneous electric field between the fermion 
bunches can be represented in terms of a linear rising po- 
tential. For larger values of the initial field strength Eq 
or the coupling e, we expect that secondary particle cre- 
ation due to the self-consistent electric field takes place. 
This mechanism would result in the depletion of the elec- 
tric field reminiscent to the effect of string-breaking. This 
will involve further studies with supercritical initial field 
strengths which is beyond the scope of the present work 
and deferred to a future publication. 



IV. CONCLUSION & OUTLOOK 

We investigated fermion-anti-fermion pair production 
in 1-|-1 dimensions based on real-time lattice simulations. 
To this end, we discussed the lattice equations of mo- 
tion using the low-cost fermion algorithm to solve them. 
In order to define gauge-invariant fermionic distributions 
corresponding to charge, energy or particlc/anti-particle 
number, we derived the lattice analogue of the contin- 
uum Dirac-Heisenberg-Wigner function. In the contin- 
uum formulation, gauge invariance of these distributions 
is achieved by a Wilson line along the straight line. On 
the lattice, however, the straight line path is not unique 
due the periodicity of the lattice. We showed that corre- 
spondence with established results is achieved by replac- 
ing the straight path in the continuum by the shortest 
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path on the lattice. 

Investigating the field-strength dependence of the 
fermion-anti-fermion production rate in a static back- 
ground field we accurately reproduced the Schwinger for- 
mula. We then discussed the decay of the field due to the 
backreaction of the created fermion-anti-fermion pairs. 
For the case of inhomogeneous gauge fields we computed 
for the first time the full problem taking backreaction 
into account. Most notably, we could show that the two 
bunches consisting of particles and anti-particles create a 
homogeneous electric field between them whereas there is 
no field outside them. In subsequent work we will extend 
these studies to supercritical initial field strengths, which 
is expected to lead to striking pair creation phenomena 
reminiscent of string breaking. 

It should be emphasized that the real-time lattice 
simulations are considerably cheaper from a computa- 
tional point of view than continuum approaches such as 
based on the Dirac-Heisenberg-Wigner function. In view 
of potential experimental applications it is crucial that 
strongly inhomogeneous configurations can be well de- 
scribed. Strong inhomogcneities are a challenge for al- 
ternative approaches based on derivative expansions un- 
derlying effective kinetic descriptions. Here the lattice 
approach, which is based on ensemble techniques using 
inhomogeneous configurations, is particularly powerful. 

We employed a low-cost fermion algorithm in our 1 + 
1 dimensional simulations even though a mode-function 
expansion of the spinors would have been the more direct 
way. One reason for our choice was that we are aiming 
at investigations of QED in 3 -I- 1 dimensions since then 
the application of the mode function expansion becomes 
impracticable. The anticipated investigations of QED 
in 3 -|- 1 dimensions will show several major differences 
compared to the massive Schwinger model. Most notably, 
the gauge degrees of freedom are dynamical in contrast to 
1-1-1 dimensions where the dynamics of the electric field is 
governed only by the fcrmionic backreaction. Moreover, 
the issue of renormalization will become relevant in 3 -f- 1 
dimensions in contrast to the super-renormalizibility of 
the massive Schwinger model. 

By means of the present lattice approach to pair pro- 
duction we hope to deepen our understanding of non- 
equilibrium QED, including highly topical issues such 
as threshold lowering, collective phenomena, instabilities 
or cascades. In the long run, this approach may also 
be extended to QCD. Most notably, the investigation 
of quark-anti-quark production from gluon fields would 
give important insights into the early stages of relativistic 
heavy-ion collisions. 
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Appendix A: Marginal distributions 

The definition of the lattice Wigner function ([62|) is 
such that we reproduce the continuum expressions for its 
marginal distributions: 

Wx{x,t) = -F{x,x;t) , (Ala) 

Regarding the position space marginal distribution on 
the lattice, we consider: 

= ^ E ■ (A2) 

Upon performing the summation over q, we encounter: 
e"^('-2')«/^= - 2N,62k,i ■ (A3) 

The Kronecker delta indicates that yV(^j) is only non- 
vanishing for even /: 

Wxil,j) = -Fil,l;j) , (A4) 

with I e A. 

Regarding the momentum space marginal distribution 
on the lattice, we consider: 

ieAw 

Due to the fact that we used the module operation in 
([S^ . we obtain: 

2Ns-l 

e-''/^=/'(fc,[/-fc]^jj) = 

(1 + e-') e-''''/^=F(fc, [I ~ fc]^,; j) . (A6) 
z=o 

The factor (1 + e""^') shows that yV(<7,i) is only non- 
vanishing for even q. Accordingly, if we redefine the sum- 
mation indices: 

li = l e A and h = [I - fclw, e A , (A7) 

we reproduce the analogue of (jAlbll : 

Wriq.j) = E E e-^"«('-'^)/^'F(?i,?2;j) , 

(A8) 

with q Cz A. 
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Appendix B: Analytic results for E{x,t) — Eq 



introduced in Sec. Ill CI are then given by: 



We briefly review some analytic results for the static 
background field [31 ■ As a matter of fact, the Dirac 
equation is analytically solvable for E{x, t) = Eq in terms 
of parabolic cylinder functions Du{z). Accordingly, it is 
also possible to compute W{x,p,t) explicitly. 

The pseudo-distributions m(x,p,t), which have been 



g{x,p,t) = , 
e{x,p,t) = [2f{p) 
'^{x,p,t) = f{p) . 



(Bla) 
(Bib) 
(Blc) 



The function /(p) is usually denotes as the single-particle 
momentum distribution: 



fip) 



,-7r/4e 



2e 



r 



-DM 



(B2) 



with 



-Diip) 

V2{P) 



|-D-l + i/2e(p)| , 
|A/2£(p)|^ , 



and 



2p 



(B3a) 
(B3b) 
(B3c) 



(B4) 



We note that /(p) is independent of the time variable t. It 
can be shown that f{p) vanishes for small momenta and 
approaches a non-vanishing constant for large momenta: 



lim f{p) = 



lim f{p) 

p— >C30 



exp I 



(B5a) 
(B5b) 



As the expressions jBTj) are spatially homogeneous, 
they are trivially related to the momentum space 
marginal distributions m-p {p, t) : 



m{x,p, t) = 



m-p{p,t) 



(B6) 



in the infinite volume L — > cx3. Most notably, the rate 
at which particles and anti-particles are created is a con- 
stant, so that the total number of particles and anti- 
particles, respectively, which are created per volume L 
and time T is given by: 



LT 



eEp 
27r 



exp 



eEn 



2tt 



exp(-^) . (B7) 
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